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We report measurements of CP violation parameters in B° — > K°ty° decays based on a data 
sample of 657 x 10 6 Bi3 pairs collected with the Belle detector at the KEKB e + e~ asymmetric- 
energy collider. We use B° — ■> K%-k° decays for both mixing-induced and direct CP violating 
asymmetry measurements and B° —+ K^ty decays for the direct CP violation measurement. The 
CP violation parameters obtained are sin20i ff = +0.67 ± 0.31(stat) ± 0.08(syst) and A K o n a = 
+0.14 ± 0.13(stat) ± 0.06(syst). The branching fraction of B° — > K°-k° decay is measured to be 
B(B° -> K°tv°) = (8.7 + 0.5(stat.) ± 0.6(syst.)) x 10 -6 . The observed A K o w o value differs by 1.9 
standard deviations from the value expected from an isospin sum rule. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw 



Decays of B mesons mediated by b — > s penguin am- 
plitudes play an important role in both measuring the 
Standard Model (SM) parameters and in probing new 
physics. In the SM, CP violation arises from a single 
irreducible Kobayashi-Maskawa (KM) phase in the 
weak-interaction quark-mixing matrix. In the decay se- 
quences T(4S*) — ► B a B° — > /cp/tag, where one of the B 
mesons decays at time tcp to a CP eigenstate fcp and 
the other decays at time itag to a final state /t ag that 
distinguishes between B° and B°, the decay rate has a 
time dependence given by 



V{At) = 



\At\/r 



4r 



B" 



1 + q- [S f sin(Am rf Ai) 
+ A f cos(Am d At)] 



(1) 



Here, Sf and Af are parameters that describe mixing- 
induced and direct CP violation, respectively, t b o is the 
B° lifetime, Arrid is the mass difference between the 
two B° mass eigenstates, At ~ tcp ~ £tag, and the b- 
fiavor charge, q = +1(— 1) when the tagged B meson is a 
B°(B°). The SM predicts S f = -£/ sin20i and A f ~ 
to a good approximation for most of the decays that pro- 
ceed via b — > sqq {q — c,s,d,u) quark transitions 
where £f = +1(— 1) corresponds to CP-even (-odd) final 
states and (f>\ is an angle of the unitarity triangle. The 
final state -^5^° is a CP eigenstate with CP eigenvalue 
£/ = —1 while K^tt is a CP eigenstate with £/ = +1. 

However, even within the SM, in B a -> K°ir° decay 
modes, both Sko^o and -4^°^° could be shifted due to 



the contribution of a color-suppressed tree diagram that 
has a V u b coupling [3| . The resulting effective parameter 
sin2<^ ff can be evaluated in the 1/rrib expansion and/or 
using SU(3) flavor symmetry 0], whereas the shift in 
Ako-jt is predicted by applying an isospin sum rule to 
the recent measurements of B meson decays into Kir fi- 
nal states [5|- The sum rule for the decay rates gives 
the following relation to within a few percent precision 
determined by SU(2) flavor symmetry [6(, 



A 



K+TT° 



2B(K+tt°)t b o 
B(K+tt-)t b+ 



B(K°tt+)t b o _ 
B(K+it-)t b+ ~ 
2B(K°tt°) 

Ak°tt° 



B(K+tt- 



(2) 



Here, B represents the branching fraction of a decay 
mode. Since the branching fractions and CP asymme- 
tries of other B — > Ktt decay modes have been mea- 
sured with good precision 0, B|, A-k jt° is constrained in 
this framework with a small error. Therefore, a signif- 
icant discrepancy between the measured and expected 
values of Ak"-k would indicate a new physics contri- 
bution to the sum rule. The expected uncertainty in 
Ako„q can be reduced by improved measurement of the 
B° — > K°tt° branching fraction. Furthermore, recent 
measurements that show an unexpectedly large difference 
between Ak+k- an d Ak+^o 0, Q makes an improved 
measurement of Ak^-k particularly interesting. In this 
paper, in addition to the B° — > K^ir mode, we measure 
the CP asymmetry in B° — > K\t^ decay for the first 
time, in order to maximize sensitivity to the direct CP 



3 



violation parameter, Ak°tt°- 

At the KEKB asymmetric-energy e + e~ (3.5 GeV on 
8 GeV) collider [l(|, the T(45) is produced with a 
Lorentz boost of /?7 = 0.425 nearly along the direction 
opposite to the positron beam line (z-axis). Since B and 
B° mesons are approximately at rest in the T(45) center- 
of-mass system (CM), At can be determined from the 
displacement in z between the fcp and /t ag decay ver- 
tices: At ~ (z C p - z tag )/(/37c) = Az/(/3 7 c). For K s tt° 
decays, the vertex position of the CP side is determined 
from the Kg decay products and the Kg mesons are re- 
quired to decay within the silicon vertex detector (SVD) 
for the time dependent CP violation measurement. Since 
we cannot obtain vertex information from K^ir decays, 
only the parameter Ak°k° is measured by comparing the 
decay rates of B° 
set of B° 

At from the decay vertex reconstruction are treated sim- 
ilarly. 

Previous measurements of CP violation in B° — > K%tt° 



Kin and B° 



K° l k°. The sub- 



Kgir events for which we cannot obtain 



decay have been reported by Belle [ll| and BABAR [K 
The previous result from Belle was based on 532 x 10 6 
BB pairs. In this report, all results are based on a data 
sample that contains 657 x 10 6 BB pairs, collected with 
the Belle detector at the KEKB operating at the T(45) 
resonance. 

The Belle detector is a large-solid-angle magnetic spec- 
trometer that consists of SVD, a 50-layer central drift 
chamber (CDC), an array of aerogel threshold Chercnkov 
counters (ACC), a barrel-like arrangement of time-of- 
flight scintillation counters (TOF), and an electromag- 
netic calorimeter (ECL) comprised of CsI(Tl) crystals lo- 
cated inside a superconducting solenoid coil that provides 
a 1.5 T magnetic field. An iron flux-return located out- 
side of the coil is instrumented to detect mesons and 
to identify muons (KLM) . The detector is described in de- 
tail elsewhere [l3| . Two configurations of the inner detec- 
tors were used. A 2.0 cm -radius beampipe and a 3-layer 
silicon vertex detector were used for the first sample of 
152 x 1Q 6 BB pairs, while a 1.5 cm -radius beampipe, a 4- 
layer silicon detector and a small-cell inner drift chamber 
were used to record the remaining 505 x 10 6 BB pairs [I! ]. 

Charged particle tracks are reconstructed with the 
SVD and CDC. Photons are identified as isolated ECL 
clusters that are not matched to any charged particle 
track. We reconstruct 7r° candidates from pairs of pho- 
tons that have energies larger than the following thresh- 
olds: 50 MeV for the barrel region and 100 MeV for the 
endcap regions, where the barrel region covers the polar 
angle range 32° < 6 < 130°, and the endcap regions cover 
the forward and backward regions. The invariant mass 
of reconstructed 7r°'s are required to be in the range be- 
tween 0.115 GeV/c 2 and 0.152 GeV/c 2 . We reconstruct 
Kg candidates from pairs of oppositely charged tracks 
having invariant mass between 0.480 GeV/c 2 and 0.516 
GeV/c 2 , which corresponds to three standard deviations 



in a Gaussian fit to the signal Monte Carlo (MC) samples. 
The flight direction of each K s candidate is required to 
be consistent with the direction of its vertex displacement 
with respect to the interaction point (IP). Candidate K\ 
mesons are selected from ECL and/or KLM hit patterns 
that are not associated with any charged tracks and con- 
sistent with the presence of a shower induced by a K\ 
meson }15j . 

For reconstructed B Q — » K s ir° candidates, we iden- 
tify B meson decays using the beam-constrained mass 



= \l ( E hc!m 



(Pb M ) 2 an d the energy difference 
AE = E% M - E™ m , where E™ m is the beam energy 
in the CM, and E% M and pg M are the CM energy and 
momentum of the reconstructed B candidate, respec- 
tively. The signal candidates used for measurements of 
CP violation parameters are selected by requiring 5.27 
GeV/c 2 < M bc < 5.29 GeV/c 2 and -0.15 GeV < AE < 
0.1 GeV. For B° — » K^tt candidates, we can only mea- 
sure the flight direction of the K®, so M^ c is calculated 
by assuming the parent B° to be at rest in the CM. The 
signal is selected by requiring M^ c > 5.255 GeV/c 2 . In 
the B° — > K s ir° analysis, if there are multiple candi- 
dates, we select the candidate that has the smallest \ 2 of 
the 7T° mass-constrained fit for the daughter photons. In 
the B° — > K^tt case, the candidate having the smallest 
cos 9 cxp is chosen, where # exp is the angle between the 
measured K° L flight direction and that expected from the 
7T momentum assuming the parent B to be at rest in 
the CM frame. 

The dominant background for the signal is from con- 



tinuum e + e 



, dd, ss or cc events. To distinguish 



the spherical BB signal events from these jet-like back- 
grounds, we combine a set of variables that character- 
ize the event topology, i.e., modified- Fox- Wolfram mo- 

iIGIIi 



ments 
variable £ s i g (bk g ) 



, into a signal (background) likelihood 
and impose requirements on the like- 
lihood ratio 7£ s /b = £sig/(£sig + £bkg): Tl s /b > 0-3 and 
K s/h > 0.5 for B° -► K s n° and B° -> K l tt candidates, 
respectively. 

The 6-flavor of the accompanying B meson is deter- 
mined from inclusive properties of particles that are not 
associated with the reconstructed B° — ► K°n° decays. 
To represent the tagging information, we use two param- 
eters, the 6-flavor charge, q and its quality factor, r 
The parameter r is an event-by-event, MC determined 
flavor-tagging dilution factor that ranges from r = for 
no flavor discrimination to r = 1 for unambiguous fla- 
vor assignment. For events with r > 0.1, the wrong tag 
fractions for six r intervals, wi (I = 1-6), and their differ- 
ences between B° and B° decays, Awi, are determined 
using a high-statistics control sample of semi-leptonic and 
hadronic b — > c decays [13, [HJ. If r < 0.1, we set the 
wrong tag fraction to 0.5, in which case the accompany- 
ing B meson does not provide tagging information and 
such events are not used for the CP violation parame- 
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ter measurement. The total effective tagging efficiency 
is estimated to be 0.29 ± 0.01, where "effective" means 
a summation over the products of tagging efficiency and 
r 2 of all types of tags used. 

The vertex position for the B° — > K s tt° decay is recon- 
structed using charged pions from the Kg decay and an 
IP constraint j2^|. Each charged pion track is required 
to have more than 1(2) hit(s) on SVD r — <j) (z) strips. 
The /tag vertex is obtained with well-reconstructed tracks 
that are not assigned to the B° — > ^571"° decay. 

Figures Q] and show the distribution of the selection 
variables for B° -> K s ir° and B° -> K° L ir° candidates. 
The signal yields are obtained from multi-dimensional 
extended unbinned maximum likelihood fits to these dis- 
tributions. The Mbc, AE and TZ B /h signal shapes for 
B° — > KgTT° are modeled with three-dimensional his- 
tograms determined from MC, while the continuum back- 
ground shapes in M^ c and AE are modeled with an AR- 
GUS function [23| and a linear function, respectively, 
whose shape and normalization are floated in the fit. The 
data from a sideband region (5.20 GeV/c 2 < Mbc < 5.26 
GeV/c 2 and 0.05 < AE < 0.20 GeV) are used to de- 
termine the continuum background shape in 7t s /b- For 
B Q — > K^ir , the signal shape in M^ c is determined 
from MC samples and the continuum background shape 
is modeled with an ARGUS function. The signal shape in 
7£s/b is determined from MC simulation. The continuum 
7\L s /b shape is determined using T(4S') off- resonance data. 
The shape of each variable in the BB background com- 
ponent is modeled using MC events. The signal yield is 
extracted in each r-bin for B° — > K^tt candidates with 
r-dependent lZ s /b shapes. For B° — > K^tt , the ratio 
of BB background to signal is evaluated from MC simu- 
lated events and the BB background contribution is then 
fixed according to that of the signal in the fit. 

We perform a fit to B° — > ^5-71"° candidates using a 
signal shape with correction factors (to account for small 
differences between data and MC) obtained from B + — > 
K + tt°. The signal yield is 634 ± 34, where the error is 
statistical only. The average signal detection efficiency 
is calculated from MC to be (22.3 ± 0.1)%. We obtain 
a B° -> K q tt° branching fraction of (8.7 ± 0.5 ± 0.6) x 
10~ 6 using only B° — > K s ir° candidates, where the first 
error is statistical and the second is systematic. The 
systematic uncertainty for the B° — » K tt° branching 
fraction is estimated by varying the correction factors 
obtained from B+ -> K+ir° by ±la (+3.6/-2.4%) and 
varying histogram probability density functions (PDF's) 
bin-by-bin by ±2cr (+1.5/— 1.6%). Uncertainties in the 
number of BB pairs (1.4%), MC statistics (0.2%), K% 
(4.9%) and ir° reconstruction efficiencies (2.8%) are also 
included. 

The signal yield of B a -> K^tt is 285 ± 52, where the 
error is statistical only. We evaluate the systematic un- 
certainty for the B° — > K^ir signal yield by smearing 
the PDF shapes of Mb c , an d r used in the fit. The 



dominant contribution is from the continuum background 
shape and the total systematic error is 20%. Taking into 
account both statistical and systematic errors, the signif- 
icance of B° -> K<lir a is 3.7ct. 




200 r 



F 1 
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R s/b 

FIG. 1: M hc -AE-H s/b fit projections of B° -> K^n candi- 
dates. The open histogram with the solid curve shows the 
fit result, the filled histogram is the BB background, and 
the dashed histogram is the sum of continuum and BB back- 
grounds. Each plot requires signal enhanced conditions for 
the other variables: 5.27 GeV/c 2 < M bc < 5.29 GeV/c 2 , 
-0.15 GeV < AE < 0.1 GeV and TZ s/b > 0.7. 
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FIG. 2: B° — * K1tv° fit projections for events with good tags: 
background subtracted Mbc (left) and lZ B / b (right). The open 
solid histogram shows the fit result. The filled histogram is 
the BB background and the open dashed histogram is the 
sum of continuum and BB background. 

We determine sin2^f ff and ^4^0 ^0 for B° — > K^ir 
by performing an unbinned maximum-likelihood fit to 
the observed At distribution. The PDF expected for the 
signal distribution, V(At; sin 20° , ^l^o^o, q, w\, Awi), is 
given by Eq. ([1]) fixing Tgo and Arrid to their world aver- 
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ages 2Jj and incorporating the effect of wrong flavor as- 
signment. The distribution is convolved with the proper- 
time interval resolution function, i? s ; g (At), which takes 
into account the finite vertex resolution. The resolution 
is determined by a multi-parameter fit to the At distribu- 
tions of hig h-statistics control samples of B° — > J/ipKg 
decays [20L [2ll| . where the K$ is used for vertex recon- 
struction. We determine the following likelihood for each 
event, 



(1 - /ol) 



/ sig 7' s ig(At')i?sig(A^ - At') 



+ (1 - f sig )V hks (At')R hke {Ati - At') 
+ f iP i(At t ). 



d(At') 



(3) 



The signal probability, / s i g , depends on r and is calcu- 
lated in each region on an event-by-event basis as a func- 
tion of Mb c ,7?. s /b and, where applicable, AE from the 
shapes given in Figs. [1] and O "Pbkg is a PDF for contin- 
uum and BB backgrounds. The background PDF's are 
determined from M^ c and AE sideband data for contin- 
uum and, MC and data for BB. The term, P D \(At) is 
a broad Gaussian function that represents a small out- 
lier component with a fraction f a \ [2C 



21|. The free pa- 



rameters in the final fits are sin 2(p1 and Ak° tt« , which 
are determined by maximizing the likelihood function 
L = U i P i (Ati;sm2(j)f l , Axo^o) where the product is over 
all events. 



The B° -> A?tt and B° 



Kgir candidates that do 



not have vertex information are only used for the deter- 
mination of Ak°tv°- Since At vanishes by integration, Eq. 
(|3|) becomes simpler: 



Pi = fsi g Psi S (q) + (1 - / S ig)7>bk g (<?), 



(4) 



where 7- > bkg('Z = ±1) = 0.5 since we assign no tag in- 
formation for the continuum background meaning that 
the number of events tagged as q = +1 and q = — 1 are 
equal. Since no CP violation is expected in the back- 
ground outlier component, we include the f \ term in the 
Pbkg PDF. The signal PDF is obtained by integrating 
the time-dependent decay rate Eq. {1} from — 00 to +00: 



1 



qAxo-^o 



1 



ioAm 2 d 



(5) 



We obtain the fit results sin2^ ff = +0.67 ± 0.31 and 
Ako^o = +0.14 ± 0.13 for B° -> A°7r°. Fits to individ- 
ual modes yield sin2</>° ff = +0.67 ± 0.31 and Ako^o — 
+0.15 ±0.13 for B° -> K%n°, and A K o^o = -0.01 + 0.45 
for B° — > A£tt°, where the errors are statistical only. 
Fig. [3] shows the background subtracted At distribu- 
tions for B° and B° tags as well as the asymmetry for 
B° — > KgTT° candidates. The dominant sources of sys- 
tematic errors are summarized in Table[I] The systematic 




FIG. 3: The top plot shows the background subtracted At dis- 
tribution for B° and B° tags where the solid (broken) curve 
represents the At curve for B° (B°) in the good tag region 
0.5 < r < 1.0. The bottom plot shows the background sub- 
tracted asymmetry defined as (N^\f - Ngf)/(Ngg + N^) in 

each At bin where N^o (A^f)) is the B° (B°) signal yield 
extracted in that At bin. The solid curve shows the CP 
asymmetry result expected from the fit. 



uncertainty from wrong tag fractions, physics parame- 
ters, resolution function, background At and background 
fractions are studied by varying each parameter by its er- 
ror. A possible fit bias is examined by fitting a large num- 
ber of pseudo-experiments. The systematic uncertainty 
for the vertex reconstruction is estimated by changing the 
charged track selection criteria. The dominant effect for 
AAk°tt° comes from misalignment between the SVD and 
CDC. The tag side interference is evaluated from pseudo- 
experiments in which the effect of possible CP violation 
in B° — » /t a g decays is taken into account [25j]. As a 
cross-check, we fit the B° lifetime using the same event 
sample that is used for the B° — > AiyTr CP violation pa- 
rameter measurement and obtain tb° — 1.46 ± 0.18 ps, 
which is consistent with the PDG world average [24j |. 

In summary, we use B° — > A<j7r° decays to measure 
the branching fraction and CP violation parameters for 
B° —> K°ir° . We use B° — > K^ir decays to measure the 
direct CP violation parameter. Our results are 



B(B° -» A°7T ) 

Ak°-k° 

sin2^ ff 



(8.7 + 0.5 + 0.6) x 10~ 6 (6) 
+0.14 + 0.13 + 0.06 (7) 
+0.67 + 0.31 + 0.08, (8) 
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TABLET: Systematic uncertainties in sin 20f and A K o n o. 



search (" Evolution of Tau-lepton Physics" 



DO 111 LC 


A «in 9rh ctt 

Z_i bill 




Wrong tag fraction 


U.UU / 


n one; 
U.UUo 




0.007 


0.001 


Resolution function 


U.Uuo 


U.UU ( 


Background At shape 


0.015 


0.006 


Background fraction 


0.029 


0.022 


Possible fit bias 


0.010 


0.020 


Vertex reconstruction 


0.013 


0.022 


Tag side interference 


0.014 


0.054 


Total 


0.077 


0.064 



where the first and second errors listed are statistical 
and systematic, respectively. These results are consis- 



tent with previous measurements ll|, [12| \ the value for 
the branching fraction is the most precise single measure- 
ment to-date. We test the isospin sum rule (Eq. [5]) by 
inserting our measured values for the branching fraction 
and Ako-k°- For the other parameters we use the most 
recent world average values [26j. We find the isospin re- 
lationship to be only marginally satisfied; the level of 
disagreement is 1.9c. Specifically, the difference between 
our measurement of Ar-o^o x B(K°tt°) and that predicted 
by Eq. [2] is 1.9a. 
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